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and shift vector) is given which preserves some of the properties of the single-hole Kerr-Schild form. 
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I. INTRODUCTION 

In the numerical simulation of gravitational spacetimes 
Einstein's equations are cast as a Cauchy problem in the 
ADM 3-1-1 formalism. In this form the equations are split 
into a set of evolution and constraint equations. The lat- 
ter set, the Hamiltonian and momentum constraints, are 
elliptic partial differential equations which must be im- 
posed at the initial slice in the Cauchy evolution. There- 
after, analytically, the evolution equations preserve them. 
York, in a series of papers gave a framework for 

solving the constraint equations known as the conformal 
formalism. We take G = 1 (Newton's constant) and c = 1 
(the speed of light). The 3-metric, g^j, is assumed to be 
of the form 

9ij = (l^^gij (1) 

where jjij is the base metric; gij is then termed the phys- 
ical metric and (f> is the conformal factor. The extrin- 
sic curvature is decomposed into its trace and trace-free 
parts: 

K^J - 4 + Ig^jK (2) 

where again the hatted quantities refer to the tensors in 
the base metric and Eij is the traceless part of the ex- 
trinsic curvature. In this approach K is assumed to be 
a given scalar function. The physical extrinsic curvature 
is related to the base extrinsic curvature through the fol- 
lowing conformal transformation: 

^y = r'"4-- (3) 

The Hamiltonian and momentum constraints then can 
be written as: 

8A(j) ^R<P- -K'^cjf' + E^^ Ei.j(j)-'^ = (4) 
3 

and 



bjE'^ - \b'K = (5) 

where A is the laplacian with respect to gij , D is the co- 
variant derivative compatible with gij and since we have 
a vacuum spacetime we take the matter terms to be zero. 

York's method, as it has been applied in the black hole 
case, assumes that the 3-space is conformally flat, with 
holes, and that the expansion of the 3-space (TrK) van- 
ishes. For the case we consider, we assume the initial 
spatial domain contains holes. A method has been given 
1^,^ for specifying an essentially analytic solution to the 
momentum constraint with symmetric boundary condi- 
tions (an infinite series has to be summed) for a multiple 
black hole spacetime. The remaining difficulty is then to 
solve the Hamiltonian constraint (Q) , which is a nonlinear 
elliptic equation, for the conformal factor. A completely 
convergent method for doing this was exhibited , where 
a controllable convergent algorithm produces solutions of 
specifiable accuracy. The static single hole initial-data 
solution found in this manner is a t ^constant slice of 
the isotropic coordinate representation of Schwarzschild 
1^ . (One of the features is that the initial slice does not 
penetrate inside the horizon.) Multiple moving spinning 
black holes can be specified. 

A new alternative to the conventional method based 
on throats and conformal imaging was proposed and im- 
plemented by Brandt and Briigmann 1^] where the black 
holes are treated as punctures. The internal asymptot- 
ically flat regions are compactified to obtain a domain 
without inner boundaries leading to significant computa- 
tional simplications. 

Recent computational black hole work has turned to 
solutions of the Kerr-Schild [|lO|,0 form, because these 
analytical models have no coordinate pathology at the 
horizon and allow slicings which penetrate within the 
horizon ||l2| . The present philosophy of handling the sin- 
gularity hidden within a black hole is to compute only up 
to the horizon [|l3|,0, though such computational tech- 
niques almost certainly require a consideration of points 
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slightly within the horizon. Further, having an analytic, 
tractable example is a great asset in developing a compu- 
tational scheme. Hence we propose a scheme where we 
utilize our freedom in setting data for the base 3-metric 
and extrinsic curvature and use this analytical specifica- 
tion based on the Kerr-Schild form. 

The motivation for this work is to initialize compu- 
tational evolutions being done by the Binary Black Hole 
Grand Challenge Alliance jl^ which has shown some suc- 
cess at evolutions based on the Kerr-Schild form. The 
approach here focuses on the 3-dimensional initial value 
problem, and puts less emphasis on a representation 
which resembles the Kerr-Schild one in a 4-dimensional 



sense. A complementary approach |16| within the Al- 
liance, which does emphasize the 4-dimensional aspects 
has recently appeared. The two approaches differ sig- 
nificantly, as we discuss more completely in section VH 
below. 



II. KERR-SCHILD FORMS FOR ISOLATED 
BLACK HOLES 

The Kerr-Schild spacetime metric is given by: 



ds^ = 7]f,^dx^'dx'' + 2H{x°')l^,l^dx^'dx'' 



(6) 



where rj^i, is the usual flat space form, iJ is a scalar func- 
tion of position and time and is an (ingoing) null vector 
(null both in the background, and in the full metric), 



with p = \/ x"^ + y'^ + z^. Notice that the function H 
given by Eq. (H) is harmonic: 



V^H = 



(12) 



where this is the flat space background Laplacian. In 
the limit r — > oo or a — > we recover the Schwarzschild 
metric in Kerr-Schild form [H = M/r, = {l,ni)). 
This corresponds to the ingoing Eddington-Finkelstcin 
|]l8| form of the Schwarzschild metric. Notice that for 
the stationary Kerr (and Schwarzschild) black holes we 
choose a particular normalization of Z^:Zt = l^t^ = 1. In 
the 3-1-1 form used in many computational approaches, 
one splits the 4- metric into a 3-metric, a lapse and a 
shift. The 3-metric gives distances measured in a given 
spacelike 3-surface (at a particular instant of time). The 
inverse of the lapse function gives the ratio of the coordi- 
nate time (measured normally to the t = const 3-space) 
to the evolution of proper time in moving along that time- 
like direction. The shift vector (multiplied by the interval 
of coordinate time dt ) gives the amount that coordinate 
labels shift in going from one t ^constant 3-space to the 
next. 

The general Kerr-Schild metric can be cast in a 3-1-1 
form: 



(lapse) a = 



1 



(shift) (3^ = 2Hlth 
(3-metric) gtj = Sij + 2Hlilj. 



(13) 

(14) 
(15) 



v'^'^U. - g'^'lJ. = 0, 



(7) 



so that If — UU. It is the fact that the Kerr-Schild form 
is based on and selects null surfaces near the black hole 
which makes it so appealing as a scheme for setting coor- 
dinates. The general Kerr-Schild black hole metric (writ- 
ten in Kerr's original rectangular coordinates) has 



H 



Mr 



a? cos^ 9 



and 



1, 



rx + ay ry — ax z 



(8) 



(9) 



Here M is the mass of the Kerr black hole and a — J/AI 
is the specific angular momentum of the black hole, r, 9 
(and (f)) are auxiliary spheroidal coordinates, z = r cos 9 
and (p is the axial angle. r{x,y,z) is obtained from the 
relation. 



x"^ + y'^ ^ z"^ ^ 
+ a^ 



(10) 



For black hole spacetimes the relation between the lapse 
and the shift obtained from this splitting guarantees that 
the horizon stays at a constant location, even though 
the lapse is non-zero at the horizon. This is in con- 
trast to the isotropic coordinate representation of the 
Schwarzschild solution, in which the shift vanishes ev- 
erywhere and the lapse vanishes at the horizon. Hence, 
for example, the static Schwarzschild black hole (the 
"Eddington-Finkelstein jl8| form") : 



H = M/r 
1 

1 + 2M/r 



2M Xi 
r r 



2M Xi Xj 



2M 



1 



v/1 + 2M/r 



r ?7jj - (2 + —)xiXj 



(16) 
(17) 

(18) 
(19) 



is smooth at the horizon r = 2M . Notice that the non- 
moving Eddington-Finkelstein metric uses 



as 



k = d^{M/H), 
where di is the partial derivative. 



(20) 



r' = -(P -a' 
2^^ 



-{p^ — a^) + a^z^. 



(11) 
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III. BOOSTED BLACK HOLES 

The Kerr-Schild metric is form-invariant under a boost 
due to its structure. This makes it ideal as a metric for 
constructing a model problem for moving black holes . 
One simply applies a constant Lorentz (boost velocity 
V as specified in the background Minkowski spacetime) 
transformation Ao to the 4-metric. The resulting metric 
retains the Kerr-Schild form, but with straightforwardly 
transformed H and I,,'. 



r'P - 



9uu 



(21) 



Under this boost, It is no longer unity. Note that because 
we start with a stationary solution the only time depen- 
dence is in the motion of the center. In this sense the 
solution has no "extra" radiation loaded into the initial 
data for the boosted hole (since we see none escape to 
infinity) . 



Under such a boost, H no longer satisfies Eq. (12), 
but every solution of Eq. ( [T^ ) in the nonmoving frame 
satisfies 



V2i/-(v-V)^H = 



(22) 



in the boosted frame, where V is the flat background 
spatial derivative operator. 



IV. BOOSTED SCHWARZSCHILD 

For example, for the Eddington-Finkelstein system 
boosted in the z-direction {v = v^), with new coordi- 
nates {x, y, z, t) we have 



r — 
k = 

— 

ly — 

Iz = 

H = 



x'^ + j^{z - vt)'^ 
7(1 — vj{z — vt)/r) 
x/r 
y/r 

7(7(2 — vt)/r — v) 
M/r 



(23) 



where 7 — 1/^ (1 — v"^) 

Under a boost, the metric becomes explicitly time de- 
pendent (because r is time dependent). 

Because the boost of the Schwarzschild solution merely 
"tilts the time axis," we can consider all of the boosted 
3-1-1 properties at an instant t = 0, in the frame which 
sees the hole moving. Subsequent times t simply offset 
the solution by an amount vt. 

After the boost, li no longer solves Eq. (EOh, but 



With Eq 
and Eq. 



I, = d^{M/H) - -fv^ . (Boosted E-F) (24) 

1^ and Eq. (M) a, are defined via Eq. (|l|) 
O). Figure (1) shows a slice through an 



Eddington-Finkelstein solution boosted to 0.5c (moving 
upward in the figure). The heavy inner contour is the 
horizon, the cardioid contours are lines of constant a, 
and the line segments indicate the direction and magni- 
tude of Pi. 

The extrinsic curvature for the boosted Kerr-Schild 
metric is given by computing gij in the boosted frame 
(given by combining the boosted terms from Eq. ( [2l| ) 
into Eqs. (|ll)-(|ll)).^ 

Now, one of the Einstein equations is 



(25) 



where the dot indicates partial time derivatives. Eq. (|2£ 
enables the computation of Kij . 
We find that 



(26) 



X \ \ \ I0.09 I / ^ 




FIG. 1. Contours of lapse and shift vector /3i for a boosted 
Schwarzschild (Eddington-Finkelstein) black hole. This is a 
cut through the equator; the hole is boosted at 0.5 upward in 
the figure. The shift is smaller ahead of the hole than behind, 
which allows the black hole to move through the computa- 
tional grid. The ovoid figure is the horizon (distorted in these 
coordinates). The cardioid curves are contours of constant 
lapse a, with values from the topmost contour down of 0.91, 
0.86, 0.78, 0.71, 0.64, 0.58, 0.52. Along the axis of motion the 
lapse for this = 0.5 case is ^3/2 at the leading point and 
1/2 at the trailing point. 



V. SETTING MULTIPLE BLACK HOLE DATA 

We will work by setting data for two Kerr-Schild black 
holes of comparable mass M\ ~ Mi- For purposes of 
development here (e.g. to estimate the size of terms), 
we also assume Afi/ri2 ^ 1 where r\2 is the coordinate 
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separation between holes. Each hole has a velocity Vi or 
V2, as appropriate, assigned to it. 
Define 

ri ^ = (x - xiYix - xiYdtj , 

r2^ ^{x-X2yix~X2y6,, (27) 

where x\ and are the coordinate positions of the holes 
on the initial slice. 

Here we use our freedom to set data and fix the back- 
ground 3-metric as follows: 

(3)c;g2 _ g^jdx^dx^ = Sijdx^dx^ + il{ri)iliiljdx^dx^ 

+ 2H{r2)2k2hdx'dx^ (28) 

Here aH , and a^i(a — 1, 2) are the functions defined from 
each single (perhaps boosted) black hole. This back- 
ground metric has two vectors corresponding to the null 
vector of the Kerr-Schild form (although at this point we 
see only their spatial components). 

The " symbol indicates that this is a conformally re- 
lated metric, while the physical metric is 

9i] = (29) 

Here (/> is a strictly positive conformal factor which will 
be determined in the process of solving the constraints. 
We start the constraint-solution process with a trial 

qKJ =ka\l)+ka\2) (30) 

These Ka ^(l)i Ka ^{2) are the individual extrinsic curva- 
tures, computed and indices raised using the single-hole 
boosted Kerr-Schild metric appropriate to either Mi or 
M2. (Henceforth we use the 2-hole physical or confor- 
mal metric.) The leading subscript on the left indicates 
that this is a zeroth order approximation (in the sense of 
"zeroth" guess; this is not an iteration method). 

Following York, we separate the trace: K = K{1) + 
K{2) from the traceless part of oKab- 

oE^ = okS - Is^.K. (31) 

K is considered a given scalar function and is not confor- 
mally scaled. The conformal scaling for qE"''' is chosen 
as: 

oE'^f' = 0-l"o^'^^ (32) 

Here qS"'' is the traceless part of the extrinsic curvature 
in the physical space associated with the zeroth guess. 
We attempt to write the momentum constraint: 

D.oEc' - \d,K ^ 0, (33) 

where Df, is the covariant derivative compatible with the 
3-metric. The momentum constraint is violated because 



of the appearance of connections from hole (1) multiply- 
ing an extrinsic curvature computed from hole (2), and 
vice versa. 

We can solve the momentum constraint equation by 
adding a term which contributes to the longitudinal part 
of the solution: 

where is a vector to be solved for 0| , and 

{Iwf^ = D^w^ + D^w'' - ^g^'^Ddw'^. (35) 

We then demand 

DbA'* - ^D^K^O, (36a) 
or 

Dbilwy'' = \d''K- DboE^''. (36b) 

This is an elliptic equation for w", which gives an addi- 
tion to the extrinsic curvature guaranteeing the solution 
of the momentum constraint. As it stands, however, we 
cannot directly solve Eq. (29) because it involves the full 
physical metric, which we have not yet specified. 
Using the result due to York Q : 

DbA''^ ^ (j)-^"^ DbA"'' , (37) 
the momentum constraint can be written as: 

mwY" = '^g^^'fbkK - DbioE-') (38) 

Here, note DbK is just dK/dx^ . 

Form (|3^) allows solution for w°' in the conformal 
frame, except that the value of (j) is not known; in fact 
this appearance of (j) leads to a coupling between this 
momentum constraint, and the Hamiltonian constraint 
to which we now turn. 

The Hamiltonian constraint solution follows York's de- 
velopment exactly. Since R — R4)^^ — 8(/)^^A(/) we have 

O = 8A0-i?0- '^K'^cjf' + (j)'\aE'^ + {IwY^) 

ioE,, + (iwYj) (39) 

Solution of the coupled set Eq. (|3|) and Eq. (|9|) for 
and w constitutes a solution of the constraint equations. 
For boundary conditions on (j) we impose ch— 1 at infin- 
ity, perhaps through a Robin condition jl^J^ , and at the 
surface horizons of the black holes (located in the back- 
ground metric) we set (j) — 1. We impose the condition 
— > at the surfaces of the black holes and at r ^ cx). 
First note that if we apply these conditions for a sin- 
gle black hole, since those exact solutions already satisfy 
the constraints we obtain a solution (j) = 1 and w'^ ~ 
immediately. Then note that if the initial configuration 
has holes widely separated, we expect ^ ~ 1 everywhere. 
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Then the momentum equation is an elliptic operation 



for w"" with source 



at large distances (since it 



arises from the Christoffel symbol x extrinsic curvature 
cross terms). Hence the "total charge" of the source is 



well localized, and w"' 



leading to corrections to 



Kab, which are the same order in r as the background 
Kab ~ r^'^. Further, for well-separated holes these cor- 
rections are small. In the Hamiltonian constraint (39), 



the extrinsic curvature terms are all squared, i.e. ~ r~ . 
Hence (again assuming ~ 1) the elliptic equation for 
(j) has finite inhomogeneous source ~ Kij ^ ^ r^^. The 
remaining question is the behavior of the Rcj) piece. How- 
ever, in the single-hole Eddington-Finkelstein metric, and 
hence in our conformal space, the behavior of R is such: 
R > 0, i? ~ at infinity, with deviations at infin- 
ity from the Eddington-Finkelstein value that go as r"**. 
Hence there is a finite contribution to the "charge" for the 
conformal factor equation arising from R, and (f> differs 
only finitely from from the Eddington-Finkelstein value 
6=1. 



VI. COORDINATE CONDITIONS FOR 
MULTIPLE BLACK HOLES 

This section deals only with the case of boosted 
Eddington-Finkelstein (non-spinning) holes. This de- 
velopment is prepared for immediate computational im- 
plementation of a multiple hole spacetime. As such, 
these ideas are tentative and need experimental (com- 
putational) verification. We present two formulations of 
coordinate setting for multiple black hole systems. Both 
are based on the idea that near the holes the space- 
time should look as closely as possible like a single hole 
Eddington-Finkelstein solution. In both cases, we set 
boundary conditions or lapse and shift near the holes 
based on the single hole analysis above. 

A defect of this current presentation is that some fea- 
tures of these coordinate specifications may be inappro- 
priate for situations with net angular momentum (e.g. 
spiraling merger that settles to a Kerr black hole), but 
these coordinate specifications will suffice for evolution a 
finite time into the future, because they are most accu- 
rate near the black holes, where the timescales associated 
with curvatures are shortest. 

We first present a scheme which is very closely linked 
to the single hole factor H. First consider the single 
boosted hole case Eq. (^, with appropriate boundary 
conditions: H = 0.5 at the horizon. H satisfies a mixed 
boundary condition at infinity. Specifying H (solving 
Eq. ( ^) for H) provides a complete solution for the gauge 
conditions, and guarantees that the algebraic relation be- 
tween lapse and shift are consistent at the horizon with 
the Eddington-Finkelstein form. Hence the horizon will 
be locked at a finite (moving, but not expanding or con- 
tracting ) coordinate location. 

Implementing this idea even in the single black hole 



case requires a kind of horizon-locking algorithm. The 
logic is: 

- locate the coordinate position of the apparent hori- 
zon. 



solve Eq. 
and a. 



2|) for H, and use H to determine Pa 



- adjust Pa to move the horizon normally back to its 
desired position. 

Iteration over these steps may be required to correctly 
generate H and Pa- 

The multiple-hole case is complicated by the fact that 

V appears in Eq. (^), and there is no uniquely defined 

V in a multiple-hole scenario. We take the point of view 
that a) the velocity of a black hole against its surround- 
ings is available by examining its recent history, b) for 
multiple hole systems the total momentum in the com- 
putational frame will be zero; where we compute the to- 
tal momentum in a completely naive way as ^ MiVj,) , 

1=1,2 

where Mi and V(j) are the parameters used in the initial 
data setting. 

We are not claiming that this will give a final con- 
figuration which is on average completely non-moving, 
but it seems likely to remove the majority of the resid- 
ual average motion. Furthermore, computational science 
is experimental science, and the results of a run show- 
ing a final drift are a signal to adjust the initial total 
momentum-setting. To apply coordinate ideas as in the 
preceding section, we would then expect that v in the co- 
ordinate specification Eq. ( p2| ) is a function of position 
having a local value appropriate to black hole 1 when at 
the horizon of black hole 1; and smoothly changing to 
the value of v associated with hole 2. A simple scheme 
would be to take the coordinate specification 



Vl 



rhi 



V2 



rh2 
r2 



(40) 



Thus velocities approach the appropriate velocity as the 
surface of the black hole is approached. The quanti- 
ties ri , r2 are computed as in Eq. ( |27| ) . The quantities 
fhi , Th2 are the horizon radii along a ray to the point we 
wish to evaluate v. The exponent p will be chosen by ex- 
periment, p = 1 gives the situation simplest to analyze. 
Then if we assume slowly moving black holes: 



2Afivi 2M2V2 
V = 1 . 

By our "zero momentum" assumption, this is 



V = 2MiVi 



1 



1 

r2 



(41) 



(42) 



For ri = r2, the point equidistant from the holes, one 
obtains v = 0. At infinity one has v ^ a.t least as r~^. 
At the surface of hole 1: 
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Vl 



2Mi 
r2 



(43) 



which gives a first order fractional correction to v. This 
may require the use of a higher positive power p in prac- 
tical applications. 

At intermediate locations we would expect H to recog- 
nize the orbital angular momentum, and take on a Kerr- 
like form, but our present development supports only the 
Schwarzschild-like case. Note, however, that the Kerr- 
like effects are "mild" until a final black hole forms, at 
which point the problem will become essential. 

We mention here a second approach which can be con- 
sidered: maximal slicing ||2^ (which determines the lapse 
a ) with minimal shear specification of the shift f3i. pl| ] 
These elliptic equations require boundary conditions on 
lapse and shift. The boundary conditions specified in 
this section determine the horizon values of the lapse and 
shift, even when we solve them via methods other than 
solving for H as suggested above. Hence those boundary 
conditions can be applied directly to the lapse and shift 
to be solved in the elliptic maximal slicing/minimal shear 
system. 



VII. DISCUSSION 

The description in section V of an algorithm to set bi- 
nary black hole data stars proceeds by generalizing the 
spatial metric to two centers, and erecting spatial "unit" 
vectors from those centers. With this form, a variant 
of standard 3-1-1 data setting is used to set up the data. 
The exposition in this paper concentrated on understand- 
ing the initial data structure of a single black hole, and 
describing how it can be used for a computational evo- 
lution. This method is to a large extent complementary 
to the other Alliance paper on setting data via the Kerr- 
Schild approach assumes a general metric (even for 
two- hole data) of the form of Eq. (^), up to terms 0{t^) 
where t = is the time of the initial spatial slice. Fur- 
thermore, the initial value problem is solved in terms of 
conditions on the null vector l^,, and on the multiplying 
scalar H, rather than expressing it directly in terms of 
the usual 3-1-1 objects, gij and Kij. While the usual Kerr- 
Schild la is shear-free, that condition is not imposed in 
this more general situation. Because of the close connec- 
tion of the approach in to the structure of the known 
Kerr-Schild solutions, a large number of cross-checks and 
analytical simplifications apply even in this more general 
situation. In part because if those simplifications, |l^ is 
able to present data for perturbed Schwarzschild black 
holes. An example set of data for "close" black holes is 
given in a perturbative limit where the deviations from 
sphericity are small. In this case then analysis can be 
carried through completely. 

Attacking the binary black hole problem computation- 
aly has sharpened our perspective that all aspects of the 



relativistic problem must be understood, for the compu- 
tation to proceed. These two approaches, or some sub- 
sequent combination, will be invaluable in setting and 
evolving binary black hole data. 
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